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Many-body states whose wave-function admits a representation in terms of a uniform binary- 
tree tensor decomposition are shown to obey to power-law two-body correlations functions. Any 
such state can be associated with the ground state of a translational invariant Hamiltonian which, 
depending on the dimension of the systems sites, involve at most couphngs between third-neighboring 
sites. Under general conditions it is shown that they describe unfrustrated systems which admit an 
exponentially large degeneracy of the ground state. 
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I. INTRODUCTION 

The selection of suitable tailored variational wavefunc- 
tions is a fundamental problem in the study of quantum 
many-body systems jL . The variational ansatz must sat- 
isfy two basic requirements: it should provide an accu- 
rate approximation of the target state (e.g. the ground 
state), and it should allow for an efficient evaluation of 
the relevant physical quantities (e.g. local observables 
and associated correlation functions). Matrix Product 
States (MPS) are a successful example of this kind 2 . It 
is possible to quantify their accuracy to approximate the 
exact wave-function j2j and in some specific cases |1| the 
ground state itself is in a matrix product form (e.g. see 
Ref. [5] for a review) . MPS are specifically suited to deal 
with not critical, short-range, one-dimensional Hamilto- 
nians. In order to overcome these limitations several gen- 
eralizations have been proposed [MS]- Projected entan- 
gled pair states [5] were introduced to deal with higher 
dimensions, weighted graph states [7] to treat systems 
with long-range interactions, and Multi-scale Entangle- 
ment Renormalization Ansatz (MERA) 8J to efficiently 
address critical systems. 

In this work we focus on one-dimensional quantum 
critical systems using homogeneous Binary- Tree States 
(HBTSs) as variational states. They share some struc- 
tural properties of scale-invariant MERA states (includ- 
ing the possibility of constructing efficient optimizing al- 
gorithms [Hllinilll]) but admit a simpler description and 
provide a prototypical realization of a real-space renor- 
malization process. Even though on general grounds it 
can be argued that these states are suitable candidates 
to approximate critical systems (e.g. they violate area 
law |12| with logarithmic corrections [13] ) an explicit 
derivation of their critical properties is still missing. We 
will prove that HBTSs can describe critical systems by 
computing the correlation functions and show that they 
decay in a power-law fashion. 

Once ascertained that HBTSs do describe critical 
ground state it would be important to know if there are 
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Figure 1: (Color online) BT network for 16 = 2" sites. Inset 
A) shows the isometric property of A; B) the maps T>l (left) 
and Vr; and C) the map S of Eq. (|3|. 



cases in which they are actually the exact ground state 
of a given model Hamiltonian. Despite the large body of 
work devoted so far on the subject there is no definite an- 
swer for critical systems (up to now only approximated, 
numerical evidences have been gathered on this issue). 
Given their ubiquitous presence of in condensed matter 
and statistical mechanics, this question is of particular 
relevance both conceptually and for possible numerical 
implementations. In this work we show that HBTSs can 
be associated, in the thermodynamic limit, to (non triv- 
ial) local and translationally-invariant parent Hamiltoni- 
ans |14j . Furthermore, similarly to what was done for 
MPS [H [TS] , we discuss sufficient conditions under which 
such operators continue to be parent Hamiltonians also 
for finite sites. By construction this allows us to iden- 
tify a class of (non trivial) unfrustrated Hamiltonians 
whose characterization has been attracting some interest 
recently, e.g. see Ref. [TBI. Although we concentrate only 
on binary trees the method we present to construct the 
parent Hamiltonian can be apply to other tensor struc- 
tures, as the MERA, which support scale invariance. 

The paper is organized as follows: Sec.[ll]is devoted to 
introduce the basic properties of HTBSs; Sec. |II B| focuses 
on how correlation functions can be computed for such 
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systems and shows the critical characters of such quanti- 



ties; Sec. Ill discuss how to construct parent Hamiltoni- 
ans for HTBSs. Conclusions and remarks are presented 
in Sec. |IV[ In the appendix we discuss how to construct 
the parent Hamiltonian in the case of a MERA state. 



II. HOMOGENEOUS BINARY TREE STATES 

Consider a ID lattice of iV = 2" sites, of a given lo- 
cal dimension d, with periodic boundary conditions. A 
generic pure state of such system can always be expressed 
as 
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(1) 



with a canonical basis for the single qudit and 

T a type-(^^) tensor. HBTSs of depth n are identified 

as those jV'^"-') whose T can be decomposed in terms of 
smaller tensors as in Fig. [l] Following Ref . |5] , each node 
of such graph represents a tensor (the emerging legs of 
the node being its indices), while a link connecting any 
two nodes represents contraction of the corresponding in- 
dices. The (yellow) element on the top of Fig. [l] describes 
a type-(^°) tensor C of elements Ci^^i^, while the 2N — 1 
triangles represent the same d x tensor A of type- (^2) 
whose elements » , satisfy the isometric condition 
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where 5" is the Kronecker delta and A"^'"^ = (A^^ ^^^)* 
is the adjoint of the A obtained by exchanging its lower 
and upper indexes and taking the complex conjugate. To- 
gether with the condition Y^e^j^ '^^1/2^'^^''^^ = 1, Eq. ^ 
automatically guarantees normalization of the HBT state 
at every level. It has been shown in [131 E] that under 
these assumptions, tensor tree states allow for an effi- 
cient evaluation of observables and correlation functions. 
In the generic case each tensor can be chosen to be differ- 
ent from the others. Being interested in scale-invariant 
systems, it is natural to assume all the tensors A to be 
equal. In the rest of the paper we will follow the formal- 
ism described in 118,. 



A. Evaluation of local observables 

a. Single-site observables. In the limit of large n, the 
physical properties of such states are fully determined by 
the Completely Positive Trace preserving (CPT) channel 
S of Fig. [l|:). It transforms a single site density matrix 
of elements [p]f = (C^IpICu) into a 2-sites states S{p) of 
elements 



(Cfl,?f2l'5(p)|C«l,^«2) = 
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Consider then a family T = {1-0^"^); ?i = 2,3, •••} of 
HBTSs of increasing sizes (depths) sharing the same A 
and C. The map S allows us to construct a simple recur- 
sive expression for the reduced density operator 
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which describes the averaged single site state of the n-th 
element of T (here /Jq"'' is the reduced density matrix of 
the a-th site of the system). Specifically the isometric 
property of A yields. 
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where V is the CPT map obtained by taking an equally 
weighted mixture of the partial traces of the map S as 
indicated in Fig. [ija. This can be expressed as 



V={Vl+Vr)/2 



(6) 



where T)l{-) = Tr2[5(-)] and Vr{-) = Tri[5(-)] with Tri,2 
representing partial traces with respect to the first and 
second site respectively. Equation ([S]) allows us to com- 
pute the average expectation of a single site observable 
0, for every full depth value n of the tree in terms of a 
repetitive application of the map V. Indeed indicating 
with (0Q.)*-"-' the expectation value on the a-th site of 
!■(/;(")) we can write 

^ f](e„)(") = Tr[e = Tr[e • P"-^Phat)] , (7) 
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where phat = Pi 
ments (^^IphatlC- 



(1) 



is the single site density matrix of ele- 
= EfcK,fcC£,fc-hC^_^Cfc,£]/2, and where 
2?" = V o ■ ■ ■ o 1) with "o" representing the composition 
of CPT maps. 

b. Two-site observables. An analogous procedure 
can be used to expressed averages of operators defined 
on z/ = 2 neighboring sites. All we have to do is to con- 
sider the density matrix 
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and build for this quantity a level-recursive mapping 
which is the two nearest-neighboring sites version of 
Eq. ([5]) (here represents the reduced density ma- 

trix of the sites a and a -\- 1 associate with a HBTS of 
depth n). The calculation is straightforward so we just 
write the result, 
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The above expression allows us to deal also with the case 
of observables operating on u neighboring sites. Indeed 
for V > Z li can be shown that any average density ma- 



trix p^/"^ can be expressed in terms of {p2"^'}m<n via the 
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application of a proper CPT map deriving from S. As 
an example, we write explicitely the case for v = 3 and 
4: 
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c. Thermodynamic limit. In the limit of infinitely 
many sites, from Eq. ([s]) it follows that if the average 
single site state associated with a HBTS of infinite depth 
characterized by a given isometry A is defined, then it 
must be a fixed point of the map V. Since CPT maps 
have a unique fixed point except for a subset of zero 
probability [12], the fixed point is defined amost-always. 
Similarly we can also provide an explicit formula for the 
thermodynamic limit of the 2-sites state (|9|, i.e. 
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This can be written either as a self-consistent equation or 
like a series in terms of by exploiting the identity ([9| 
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The series in convergent in any norm, thanks to the ge- 
ometric factor and the fact that CPT are non expansive. 
Such argument becomes even simpler when dealing with 
three or more n-n sites density matrices. Indeed one can 
show that for any integer z/ there exists a CPT map X'2-j-i/ 
such that, the v nearest neighbors sites density matrix 
/oi,°°^ (averaged over translations) in the thermodynamic 
limit is given by. 
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This provides a complete characterization of the local 
properties of our infinitely deep HBTS. For future refer- 
ence we report the expression for case v — 3 and 4, 

2?2^3 = {Vr<S)S + S<S)Vl)/2 , (15) 
2?2^4 = (5(g)5+(I?i?®5(g)I?L)o2?2^3)/2,(16) 

(notice that 2?2->-3 is exactly the channel which enters in 
Eq. (|lO|)). Finally we notice that all these quantities are 
independent from the element C of the HBTS, implying 
that in the thermodynamical limit, the local structure of 
the state loses all its dependence from such element. As 
the physics of the system is determined by the algebra of 
the local observables, this implies that all HBTS of infi- 
nite depth, associated with a given A but with different 
C describe the same state. 



B. Correlation functions 

Similarly to what has been done for the local observ- 
ables in the previous section, also correlation function can 



be expressed in terms of iterative application of certain 
maps. Most important for our work is to show that this 
procedure leads naturally, in the case of homogeneous 
trees, to power-law decays for the correlators: the expo- 
nents being related to the eigenvalues of the map |18j . 

In this following we focus on two-point correlation 
functions. As discussed before, since HBTSs are not man- 
ifestly translationally invariant, an average over transla- 
tions has to be made. We thus introduce the quantities 
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with 8 and O' being two single sites observables. A re- 
markable simplification is achieved for any distance Aa 
equal to a power of 2. Under this condition we find that 

ci;l2™ = Tr[(e ® e')^'"(pt""^ - v[7"'^)] , (17) 

where p is the map ;0 = ^{Vl^Vl+Vr® Vr). The 

quantity r]^ { is the averaged 2 sites nearest neighbour 
density matrix after we traced away every quantum cor- 
relation, while keeping eventual classical correlations in- 
tact, i.e. 
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Take then n ^ 00 while keeping m = log2 Aa fixed. In 

this context it is important to notice that, like P2"'' ^ilso 

rji { has a well-defined limit. It coincides with the two 
sites state, 
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with X^being the fixed point of Exploiting this fact we 
can thus write the thermodynamic limit of the correlation 

as 

= TT[{p^r^ -n^{f) X"(e®e')], (20) 

where is the adjoint superoperator of ^ (with respect 
to the operator scalar product {A,B) = Tt[A'^B]). Re- 
call that we are keeping track of the quantum correlations 
^(00) _ fj^^°^^ is a traceless matrix, 

-Am! ~^ fo'^ ~^ because 
;0°^{X) = XTr[X]. This shows that the only residual 
influence on m is kept through the number of times we 
have to apply the appropriate map to its matrix argu- 
ment. By decomposing jA in Jordan blocks, one finds its 
set of eigenoperators; let us assume that Q (E) O' is one 
of such operators, related to the eigenvalue k, then the 
correlation function is expressed as follows 
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only, moreover, since P2 

we are guaranteed that 
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where we separated the prefactor 
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(22) 



The critical exponents are defined by the spectrum of 
and the related primary fields are given by the respec- 
tive eigenoperators. Notice that such exponents have al- 
ways negative real-part, since all |k| < 1 because is a 
completely positive unital operator (and with the mix- 
ing condition only |At| < 1 and k = 1 are allowed, e.g. 
see Ref. [10]); this guarantees that such correlations are 
actually decaying power-law functons. 

If the observables 9 (g) 6' are not an eigenoperator of 
_/( their correlator is typically a sum of power-laws and 
may manifest logarythmic corrections (arising from the 
fact that jA is generally not diagonalizable) 
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where the sum spans over the spectrum of y4', and (?fc(-) 
are polynomials in their main argument. The present 
considerations prove the criticality character of HBTS. 



III. PARENT HAMILTONIANS 

In the previous section we showed that HBTSs support 
power-law decay of correlators and we related the asso- 
ciated critical exponents to the tensors which the define 
the state. Is there any case where a HBTS is the exact 
ground state of a short-range critical Hamiltonian? In 
this section we show how to construct local translation- 
ally invariant Hamiltonians for which a given homoge- 
neous HBTS is the exact ground state. First we focuses 
on the case of infinite dimensional systems (thermody- 
namical limit). Then we show how the analysis can be 
extended to the case of finite dimensional HBTSs. 



A. Thermodynamic limit 

Consider a Hamiltonian which involves at most (j^— 1)- 
neighboring couplings of the form 
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The factor 1 /N is introduced to keep a finite spectrum in 
the thermodynamical limit, and where H^(a) is an inter- 
action term that couples v consecutive sites starting from 
the a-th (i.e. the sites a, • • • jV—l + a). The expectation 
values over the infinite HBTS of this Hamiltonian is 



(25) 

with pl^^ (the averaged I'-neighboring sites density ma- 
trix) being a quantity we can calculate as discussed in 
the previous sections. 



Let us for a moment assume that the rank of /ol°°'' is 
less than its maximum d'^ . This means that such density 
matrix has a nontrivial kernel, which can be decomposed 
in a basis of vectors {\(p^{k))}k- Therefore we take 



(26) 



with Ek being arbitrary positive constants. This is pos- 
itive by construction, and so is the associated H. Then, 
since the image of H,^ belongs to the kernel of , it is 
clear that pi°°^ = 0, and so = as weU. In 

the end, we built a positive, translation invariant, Hamil- 
tonian, with {v — l)-neighboring coupling terms, whose 
expectation value over our HBTS is zero; this means that 
the state is a ground state for H. The only caveat to make 
it work is to demonstrate that, for some u we have 



rank(p(,°°)) < , 



(27) 



(otherwise would be the trivial null operator). The 
fundamental ingredient to verify this is to notice that 
the channel S of Eq. ^ preserves rank while increas- 
ing dimensions (i.e. it is an isometric mapping). Let 
thus investigate the case = 3. We know that the state 
^3°°^ is obtained by exploiting the first of the mapping of 
Eq. (16). Specifically we have 
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{I(g>S)iA) + {S(g>X){B) , (28) 



with X being the single site identity mapping and with A 
and B some d^ x d^ positive matrices. The maps I (E) S 
and S <ET preserve the rank, and the rank of the sum is 
less or equal than the sum of ranks, thus leading us to 
the inequality 



rank(/93°°^) < 



2d:' 



(29) 



over a maximum of d^. Therefore if the local dimension 
d is 3 (spin 1) or greater then we already achieved our 

goal of finding a pi°°^ matrix with non-maximal rank. 

For d — 2 (spin 1/2) instead we have to move to = 4. 
In this case the state to consider is 
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= (5®5)(A')-H(X®5(»T)(B') • (30) 



Since its rank obeys the inequality 



rank(p 



(00)n 



<d' 



(31) 



we have found a state that already for d = 2 possess 
a nontrivial kernel (indeed in this case rank(y04°°^) < 12 
which is strictly minor than the total dimension d"* = 16). 

In summary this shows that any given infinite HBTS 
admits always a local translationally invariant non-trivial 
parent Hamiltonian "H, which can be constructed explic- 



itly as in Eq. (26). For o? > 3 such Ti. can be chosen to 
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have interactions which involve up to second neighbor- 
ing couphngs. For d = 2 instead we can always choose 
T-l with up to third neighboring couplings. Our analy- 
sis deals with the worst case scenario, if it occurs that 
P2°°^ is non full rank by accident, one can construct a 
shorter-ranged (i.e. nearest neoughbours) parent hamil- 
tonian using ( 26 ) . 



B. Finite size case 
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The above approach formally applies to the case of 
infinitely many sites, and in general there is no guarantee 
that the selected H will be a parent Hamiltonian of the 
HBTS IV'^"^) when n is finite. Nonetheless the proof can 
be extend to cover also this case in most situations. This 
will allow us to prove that for N even, T-L is unfrustrated 
and that its ground space must have dimension Dgr larger 
than 

To show this we focus on the case d > 3 and assume 
that our HBTS has P2°°'' ^^1^ rank (generalization to 



d = 2 shall be dealt later) ; this guarantees that Eq. ( 26 ) 



provides a parent hamiltonian for the thermodynam- 
ical state with a three-body interaction H^. Consider 
then a generic state \ip) of N/2 sites and "grow" a BT 
level from it, using the same A isometry we used to build 
H. This way we obtain a A^-sited state 



(32) 



which, by varying spans a subspace 6 of dimension 
(when N is power of 2 one element of such sub- 
space is for instance the HBTS we started with). The 
expectation value {(j)\'H\4>) can then be expressed as 

Tr[g3 H3] = Tr[2?2^3(r2) i/3] = Tr[f2 ^3^2(i?3)] , 

where 53 is the averaged reduced density matrices of 3- 
neighboring sites of {(p), r2 is the averaged reduced den- 
sity matrices of 2-neighboring sites of \ip), while A3^2 is 
the Heisenberg conjugate map of I'2->-3- At this point we 
observe that ^3-^2 (.^3) is the null operator. This follows 
form Eq. ([T6| which allows us to write 



= Tiipr' H,] = Tr[2?2^3(p^°°^) H,i 

>2(H3)], 



Tr[pi°°U: 



(33) 



where the first identity simply states that H is the 
parent Hamiltonian of the HBTS at thermodynamical 
limit. Since p^^^ has maximal support by hypothesis 
and ^3-5.2(^3) is positive semidefinite by construction, 
Eq. ^ implies ^^3^2(^3) = 0. Equation (fSSl) then 



leads to (^iT^jc/)) = which, together with the positiv- 
ity of H, tells us that each one of the vectors {(f)) of the 
subspace 6 is a ground state of the parent Hamiltonian 

n. 

Let us now deal briefly with the case d — 2. If ^3°°"* 
is nonfuU rank then we can build a three-body interact- 
ing parent hamiltonian just like the c? > 3 case, and the 



Figure 2: (Color online) Unnormalized density of states of the 
parent Hamiltonian generated from a sample HBTS for A'' = 8 
sites (energy levels have been rescaled to the maximum energy 
eigenvalue). For this example it occurs that the ground state 
degeneracy is twice the lower bound we discuss in the paper: 
32-fold over 256 states. More precisely, we have shown that 
the ground space of the Hamiltonian in this case coincides 
with the direct sum © © T{<S) of the space & formed by the 



vectors defined in Eq. (321, and by its translated version r((S), 



T being the one-site translation (we verified that in this case 
the two sectors form linearly independent subspace). 



generalization to the finite setting is identical. Other- 
wise, via P4°°'' we can build a positive parent hamilto- 
nian Ti, of the thermodynamical state with four-body in- 
teractions 7?4. Evaluating its expectation value on the 



iV-sited state ( 32 1 it is then easy to verify that it nullify 



(the proof is similar to the previous case, and it exploits 
the fact that p^^^ has full rank). 

The above discussion proves that for all even N , the 
Hamiltonian % (resp. %') has a ground eigen-space 
which is at least d^/^ dimensional PTI . The presence of a 
wide ground state degeneracy is in accordance with sym- 
metry predictions: since a finite HBT state |'0) breaks the 
translational symmetry at every lenghtscale, the whole 
space generated by ,T\ip) ,T'^\ii)) . ..T^\ip)} must be 
embedded within the ground space. The present argu- 
ment also implies that % represents an unfrustrated sys- 
tem. Indeed if {4>\'H\(j)) = then each local component of 
Ti, needs to nullify on i.e. {(j)\H^{a)\(j)) — Va. As an 
example in Fig. [2] we report the eigenvalues degeneracies 
for a parent Hamiltonian % generated from an isometry 
A defined by the following mapping 



|0) ^ 101) , 

|1) ^ ^(100) + 111)), 



(34) 
(35) 



(here d = 2 while Ti, was generated by taking the free- 
parameters Ek of Eq. (26) to be uniform). For N — 



4, 6, 8 the ground state degeneracy turns out to be exactly 
2 2^/^ showing that in this case 6 and T(6) [5T] saturate 
completely the corresponding eigenspace (the figure only 
reports the case TV = 8). We also checked numerically 
the case of N odd (for which the previous theoretical 
analysis does not hold): in this case the ground state 
energy is not null showing that % is frustrated and that 
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its degeneracy is way smaller than d^/^, and controllable 
by choosing the Ek parameters appropriately. 

This procedure to construct a parent Hamiltonian for a 
binary tree can also be applied to other tensor structures. 
In the Appendix we discuss the case of a MERA. 



scales with the square root of the dimension of the Hilbert 
space. Similar results can be obtained (see the appendix) 
also for MERA states. 
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IV. CONCLUSIONS 

In this paper we analyzed the potential of binary trees 
to simulate efficiently quantum critical systems. Previ- 
ous works [T3] presented evidences that such states could 
yield violation of the area law with logarithmic correc- 
tions. In this paper we focused on homogenous configu- 
rations which allow for an explicit analytic treatment of 
the thermodynamical limit. Their hierarchical, scale in- 
variant structure suggest that they should be capable of 
exhibit critical behaviors, at least once a proper averag- 
ing over translations has been performed to compensate 
for their explicit lack of translational invariance [53] . For 
instance by looking at their tensorial decomposition it 
is clear that HBTS do not violate the area law for all 
possible partitions of the sites (e.g. since the left side of 
the graph is connected with the right side by only a sin- 
gle link, the resulting block entropy will be independent 
from the number of sites). It is reasonable to belief that 
such "anomalies" however will wash away when averag- 
ing over all possible translations (a legitimate operation 
when simulating translationally invariant systems). In 
the case of the block entropy this can be heuristically 
verified by noticing that indeed the average number of 
tensor links that needs to be cut in order to disconnect 
the causal cone of a block of consecutive sites from the 
rest, scales almost logarithmically with the block size. 
To test the validity of these arguments, in our paper we 
focused on the behavior of two-point correlation func- 
tions in the thermodynamical limit of infinitely many 
sites. Once averaged over all possible translations, we 
proven that these quantities can be explicitly computed 
and showed that they decay as power law in agreement 
with the criticality character of HBTS. 

In the second part of the paper we have then shown 
that HBTS are the exact ground states of short-range 
interacting Hamiltonians. In particular we gave a pro- 
cedure to built such parent Hamiltonian. In the case 
we analyzed the ground state has a degeneracy which 



We acknowledge fruitful discussions with G. E. San- 
toro and financial support from IP-EUROSQIP, FIRB- 
RBID08B3FM, SFB/TRR 21, and the National Research 
Foundation and Ministry of Education Singapore. 
Appendix A: Parent Hamiltonian for MERA states 

In this Appendix we discuss how to generalize the anal- 
ysis of Sec. HI to the case of scale-invariant (i.e. homo- 
geneous) MERA states ^8,. Indeed also for a MERA it 
is possible to establish upper bounds for the rank of the 

states (the translationally averaged neighbouring 
sites density matrix in the thermodynamical limit) by 
exploiting growth superoperators properties, for varous 
block sizes v. Therefore, by finding the suitable (small- 
est) V for which such rank is not maximal, the construc- 
tion of a parent hamiltonian interaction term according 
to equation ( 26 ) is straightforward. 



Precisely, such minimal parent interaction lenght v de- 
pends on the topology of the original MERA [22] and 
its local dimension d. For a binary MERA structure, 
we find that rank(pg°°'') < 2d^ whose value is not max- 
imal for d > 3; for completeness rank(pg°°'') < d"^ + 
takes care of the case d = 2; thus is always possible to 
build a parent hamiltonian with a 5 or 6 bodies interac- 
tion. When considering to a ternary MERA structure, 
we have to involve a block of seven nearest neighbours to 
achieve the meaningful bound rank(y07°°^) < 3d^, always 
nonmaximal regardless of d. Furthermore, one can pro- 
duce analogous conditions under which "H will be still a 
parent (unfrustrated) Hamiltonian for the finite site case 
and verify that it posses a ground state degeneracy which 
is exponentially large (order d^/^ or d^^^). A main dif- 
ference between this case and the previous one is that 
a parent hamitonian for a finite HBT state also admits 
always a dimerized ground state (just pick up a vector of 
Eq. (32) that is build by taking lip) as a product state), 
while in the MERA context there is no such proof of 
triviality. 
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